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The calculation of neutrino decoupling from nuclear matter requires a transport formalism capable of handling 
both collisions and flavor mixing. The first steps towards such a formalism are the construction of neutrino and 
antineutrino 'distribution matrices,' and a determination of the Liouville equations they satisfy in the noninter- 
acting case. These steps are accomplished through study of a Wigner-transformed 'density function,' the mean 
value of paired neutrino quantum field operators. 



A massive star supports itself against gravity 
through the energy generated by a succession of 
nuclear fusion reactions; but once the core burns 
to iron-group elements it can extract no further 
energy from fusion, and when electron degeneracy 
can no longer support it, catastrophic collapse re- 
sults. 

The halt of collapse results in a shock wave 
that will eventually disrupt the entire star. The 
inner part of the core collapses subsonically until 
it reaches nuclear density or so; but the outer part 
of the core collapses supersonically, and it doesn't 
receive word that the collapse of the inner core 
has halted until it slams into the inner core and 
forms a shock wave. The shock moves outwards 
before stalling due to energy losses from neutrino 
emission and the breakup of heavy nuclei falling 
through the shock. 

The mechanism of shock revival is still un- 
known, but simulations exploring the mechanism 
have to take account of the neutrino heating and 
cooling rates behind the shock. These rates de- 
pend on the energy and angle distributions of neu- 
trinos, which must be tracked as the neutrinos 
decouple from the newly-born neutron star. 
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Microscopic quantum mechanical effects can be 
included if the tracking of neutrino energy and an- 
gle distributions is handled with classical trans- 
port. A classical distribution function /(t,x, p) 
gives the differential number of neutrinos with 
momenta within (fip of p at positions within d^x 
of x: 



dN = f{t,x,p) ^ d'^. 



(1) 



(For neutrinos the spin degeneracy g is 1). The 
evolution of / is given by the Boltzmann equation 
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The left-hand side is the Liouville operator. It is 
closely related to the geodesic equations describ- 
ing classical particle trajectories in phase space. 
The first term gives the spacetime dependence, 
and the second term includes (for instance) par- 
ticle redshifts and trajectory bending due to gen- 
eral relativity. The collision integral on the right- 
hand side describes particle interactions that add 
or remove particles from a particular trajectory. 
Rates computed with (for instance) quantum field 
theory can be used here, so that inherently quan- 
tum mechanical — but microscopic — phenomena 
such as particle creation and annihilation, and 
phase space blocking, can be included. 
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However, neutrino flavor mixing is a macro- 
scopic quantum effect that cannot be treated with 
classical transport, which refers only to particle 
number and not amplitudes or phases. 

What about the evolution of quantum occupa- 
tion numbers rather than classical transport? An 
occupation number n{t, q) gives the average num- 
ber of particles occupying momentum eigenstate 
q in the (effectively infinite) quantization volume 
V. This is most naturally useful in the case of 
spatial homogeneity. The time derivative dn/dt 
is given by quantum transition rates. 

Noting that neutrino interaction rates have 
off-diagonal elements when written in terms of 
'physical states' (or mass eigenstates) , Dolgov 
[1] expanded the occupation number to a matrix 
p{t, |q|) in flavor/mass space in a study of neu- 
trinos in the (homogeneous and isotropic) early 
universe. 

Keep in mind the differences between classical 
distribution functions /(i,x, p) and quantum oc- 
cupation numbers n(t, q). The momentum of a 
classical particle is here labeled by p, while a mo- 
mentum eigenstate is labeled by q. In the classi- 
cal limit of a quantum description, p might be the 
centroid of a wave packet superposition of eigen- 
states q. 

However, in evolving the neutrino occupation 
matrix in the context of the early universe, the 
classical Liouville operator was employed to ac- 
count for the redshift associated with cosmic ex- 
pansion. That is, Dolgov tacitly assumed the sub- 
stitution 
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in which the time derivative of the quantum occu- 
pation matrix goes over to the Liouville operator 
acting on a distribution matrix whose momentum 
dependence has its classical meaning. 

Unlike the early universe case, neutrino trans- 
port in supernovae requires treatment of spatial 
inhomogeneity. The 'Liouville replacement' of 
Eq. (I3|) seems natural enough in the homogeneous 
case, but the introduction of spatial dependence 
should give one pause, because of the quantum 
mechanical incompatibility of position and mo- 
mentum. Nevertheless, in the limit in which the 



neutrino trajectories are classical and only the fla- 
vor evolution is quantum mechanical, one expects 
the same terms as in the classical Liouville oper- 
ator to emerge, including the spatial derivative 
term. 

The Liouville operators for neutrino and an- 
tineutrino distribution matrices also contain ef- 
fects of flavor mixing. In flat spacetime, and in 
the absence of interactions, the distribution ma- 
trices are governed by 

p'^^p(i,x,p)-f ^[A,p(t,x,p)] = 0, (4) 

p^^p(i,x,p)-^[A,p(t,x,p)] = 0. (5) 

Species indices are suppressed. Only squared 
mass differences affect observable oscillation phe- 
nomena, and so a squared mass matrix in the 
commutator has been replaced by 



A = Af2 - -Ti(M 



(6) 



The commutator enters with opposite sign in the 
case of antineutrinos if one wants the neutrino 
and antineutrino distribution matrices to trans- 
form in the same way under a flavor/mass basis 
change. 

How can the Liouville equations of Eqs. ([4]) and 

— in particular, their spatial dependence — be 
derived rather than just assumed? 

One way is through a covariant version of the 
familiar simple model of ffavor mixing [5]. In 
this version a Schrodinger equation for a neu- 
trino state is written down that involves not evo- 
lution in time or space alone, but along a world- 
line with affine parameter A. One can build a 
density operator out of the neutrino states; the 
Liouville equation then follows from application 
of the Schrodinger equation to the kets and bras 
out of which the density operator is constructed, 
together with the relation d/dX — p^^d/dx^ . 

But we will ultimately want a formalism in 
which collisions can be included, and as a first 
step towards this, the Liouville equations can be 
constructed from the equations of motion for a 
quantum 'density function.' the mean value of 
paired neutrino quantum field operators [2]. 

Define a density function r|™(?/, z) as the mean 
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value of a pair of normal-ordered neutrino quan- 
tum field operators: 

irff(y,z) = (7V^f(y)P™(z)). (7) 

It depends on spacetime positions y and z, and 
has mass indices i and j and spinor indices ^ and 
m. The A'' denotes normal ordering. Writing the 
quantum field operator as v{y) = A{y) + B{y) in 
terms of its positive- and negative- frequency parts 
A{y) and B{y), the density operator becomes 

ir^™(y,z)^-(if(z)Af(y)) + (i3f(2;)Bf(z)),(8) 

in which we have separate neutrino and an- 
tineutrino terms (rapidly oscillating cross terms 
not relevant to the macroscopic limit have been 
dropped) . 

The mixed representation of the density func- 
tion begins to look more like a classical distri- 
bution function. Introduce the average position 
variable x — (jj + z)/2 and the spatial difference 
variable S = y — z. The mixed representation 
is obtained with a Wigner transformation, which 
is a Fourier transformation with respect to the 
difference variable: 

g^rix,P)^ I d^E e'^-rff (^x-H|,a:-|).(9) 

Separate neutrino and antineutrino distribu- 
tions G'f™(x,p) and Gf™(a;,p) are obtained from 
Qf^{x,P) by projecting out its positive- and 
negative- frequency parts. 

In general circumstances the mixed represen- 
tation provides only complementary position or 
momentum distributions, when integrated over 
momentum space or the volume of the system re- 
spectively. This is in accordance with the quan- 
tum mechanical incompatibility of position and 
momentum. 

However, a joint position/momentum distribu- 
tion does emerge under classical conditions. This 
can be shown by explicitly evaluating the mixed 
representation with respect to a state of N uncor- 
related wave packets [2]. After calculation, and 
arguments about coarse graining, and application 
of classicality conditions of weak inhomogeneity 
and slow time variation, a neutrino distribution 



matrix built on classical trajectories emerges: 

N 

p,j{t, X, p) = ^(27r)3 _ x,(t)) 6Hp - Pa) 

a=l 

xe-'(^p»--^p»-)*{/„„,[/*^,. (10) 

(The antineutrino distribution matrix is the same 
except for a sign reversal in the exponential.) 
Here Xo(t) is the classical trajectory followed by 
the ath neutrino, whose wave packet centroid is 
Pa. This expression corresponds to a system of 
neutrinos born in flavors Oa at t — 0; the ener- 
gies corresponding to the superposed mass eigen- 
state wave packet centroids are Ep^^i. The trans- 
formation to a flavor-basis distribution matrix 
follows directly from Eq. ^ and the relation 
i^aiy) = Uai Viiy) between flavor and mass fields. 
The diagonal elements of the flavor-basis distribu- 
tion matrix are just as expected for a collection of 
neutrinos whose flavor content varies along classi- 
cal trajectories according to the familiar vacuum 
oscillation probabilities. 

The Liouville equations follow from the differ- 
ence of Klein-Gordon equations 

(□, - □.) r^'"(y, z) + [A, r^™(y, z)] = o (ii) 

obeyed by the density function. This can be seen 
by noting that Dy — Dz — 2 ■ in the co- 
ordinates of the mixed representation: when Eq. 
(Ilip is transformed by the inverse of Eq. ^ and 
its positive- and negative- frequency parts are pro- 
jected out, Eqs. (HI and ([5]) are the result. 

The calculations outlined here serve as a first 
step towards the handling of neutrino interac- 
tions with a diagrammatic approach based on a 
nonequilibrium Green's function [3]. In addition 
to a Green's function, this approach — which is 
sometimes called 'Keldysh theory' — involves the 
density function considered in this paper and two 
other types of field operator pairings. 
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